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INTRODUCTION 


This  annual  summary  report  describes  research  performed  from  1  August  1989 
through  30  September  1990  (14  months)  with  support  from  ONR  under  grant 
N00014-89-J-1003.  The  following  projects  are  discussed  in  this  report: 

I.  Pulsed  Finite  Amplitude  Sound  Beams 

II.  Nonlinear  Acoustic  Wave  Propagation  in  a  Liquid  Layer 

III.  Nonlinear  Effects  in  Asymmetric  Cylindrical  Sound  Beams 

IV.  Effects  of  Absorption  on  the  Interaction  of  Sound  Beams 

V.  Parametric  Receiving  Arrays 

Contributions  to  these  projects  were  made  by  the  following  individuals: 

Senior  Personnel 

•  M.  F.  Hamilton,  principal  investigator  (projects  I-V) 

•  J.  Naze  Tj0tta,  visiting  scientist  (projects  IV  and  V) 

•  S.  Tj0tta,  visiting  scientist  (project  IV) 

•  E.  A.  Zabolotskaya,  visiting  scientist  (project  II) 

Graduate  Students 

•  C.  M.  Darvennes,  Ph.D.  student  in  Mechanical  Engineering  during  August- 
December  1989,  postdoctoral  research  fellow  during  January- August  1990 
(projects  IV  and  V) 

•  E.  E.  Kim,  M.S.  student  in  Mechanical  Engineering  (project  III) 

•  Y.-S.  Lee,  Ph.D.  student  in  Mechanical  Engineering  (project  I) 

Professor  Naze  Tj0tta  spent  the  entire  14  months  covered  by  this  report,  and 
Professor  Tj0tta  spent  four  months,  at  Applied  Research  Laboratories  at  the  Uni¬ 
versity  of  Texas  at  Austin  (ARL:UT)  while  on  leave  from  the  University  of  Bergen, 
Norway.  The  Tj0ttas  received  no  support  from  the  ONR  grant  covered  by  this 
report.  Their  support  on  the  above  projects  was  provided  by  AFL:UT,  the  Texas 
Advanced  Research  Program  (TARP),  VISTA/STATOIL  of  Norway,  and  the  Odd 
Hassel  Fellowship  of  the  Norwegian  Research  Council  for  Science  and  Humanities 
(only  JNT).  Professor  Zabolotskaya  spent  five  months  in  the  Mechanical  Engineer¬ 
ing  Department  while  on  leave  from  the  Institute  of  General  Physics  in  Moscow, 
USSR.  Although  she  received  her  entire  support  from  the  David  and  Lucile  Packard 
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Fellowship  for  Science  and  Engineering,  her  research  was  performed  in  collabora¬ 
tion  with  Hamilton,  whose  work  was  supported  in  part  by  ONR.  Kim  received 
partial  support  from  TARP  (2  months),  while  Lee  received  his  primary  support 
from  TARP  (11  months).  Hamilton  also  received  partial  financial  support  from 
TARP. 

The  following  manuscripts  and  abstracts,  which  describe  work  supported  at 
least  in  part  by  ONR,  have  been  published  (or  submitted  for  publication)  since 
1  August  1989. 

Refereed  Publications 

•  C.  M.  Darvennes  and  M.  F.  Hamilton,  “Scattering  of  sound  by  sound  from 
two  Gaussian  beams,”  J.  Acoust.  Soc.  Am.  87,  pp.  1955-1964  (1990). 

•  C.  M.  Darvennes,  M.  F.  Hamilton,  J.  Naze  Tj0tta,  and  S.  Tjptta,  “Effects 
of  absorption  on  the  nonlinear  interaction  of  sound  beams,”  submitted  in 
April  1990  for  publication  in  J.  Acoust.  Soc.  Am. 

•  M.  F.  Hamilton  and  D.  T.  Blackstock,  “On  the  linearity  of  the  momen¬ 
tum  equation  for  progressive  plane  waves  of  finite  amplitude,”  submitted  in 
June  1990  for  publication  in  J.  Acoust.  Soc.  Am. 

•  M.  F.  Hamilton  and  E.  A.  Zabolotskaya,  “Nonlinear  acoustic  wave  propa¬ 
gation  in  a  liquid  layer  between  a  rigid  and  a  free  surface,”  submitted  in 
June  1990  for  publication  in  J.  Acoust.  Soc.  Am. 

Publications  in  Conference  Proceedings 

•  C.  M.  Darvennes,  M.  F.  Hamilton,  J.  Naze  Tj0tta,  and  S.  Tj0tta,  “Scattering 
of  sound  by  sound:  Effects  of  absorption,”  in  Proceedings  of  the  13th  Interna¬ 
tional  Congress  on  Acoustics,  Belgrade,  Yugoslavia,  August  1989,  edited  by 
P.  Pravica  and  G.  Drakulic  (Sava  Centar,  Belgrade,  1989),  Vol.  1,  pp.  283- 
286. 

•  M.  F.  Hamilton,  “Finite  amplitude  sound  in  waveguides,”  Compilation  of  Ab¬ 
stracts,  XIII  Scandinavian  Cooperation  Meeting  in  Acoustics/Hydrodynamics, 
edited  by  H.  Hobaek  (Scientific/Technical  Report  No.  227,  Department  of 
Physics,  University  of  Bergen,  Norway,  1990),  pp.  4-9. 

•  C.  M.  Darvennes,  M.  F.  Hamilton,  and  J.  Naze  Tjtta,  “Parametric  reception 
near  a  reflecting  surface,”  Frontiers  of  Nonlinear  Acoustics:  12th  ISNA,  edited 
by  M.  F.  Hamilton  and  D.  T.  Blackstock  (Elsevier,  London,  1990),  pp.  239 
244. 
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•  M.  F.  Hamilton  and  E.  A.  Zabolotskaya,  “Nonlinear  wave  propagation  in  a 
fluid  layer,”  Frontiers  of  Nonlinear  Acoustics:  12th  ISNA,  edited  by  M.  F. 
Hamilton  and  D.  T.  Blackstock  (Elsevier,  London,  1990),  pp.  321-326. 

Oral  Presentation  Abstracts 

•  C.  M.  Darvennes,  M.  F.  Hamilton,  J.  Naze  Tjotta,  and  S.  Tjotta,  ‘  Effects 
of  focusing  on  the  scattering  of  sound  by  sound,”  J.  Acoust.  Soc.  Am.  86, 
S106  (1989). 

•  C.  M.  Darvennes,  M.  F.  Hamilton,  and  J.  Naze  Tjptta,  “Parametric  reflection 
near  a  reflecting  surface,”  J.  Acoust.  Soc.  Am.  87,  S21  (1990). 

•  M.  F.  Hamilton  and  E.  A.  Zabolotskaya,  “Nonlinear  wave  propagation  in  a 
fluid  layer,”  J.  Acoust.  Soc.  Am.  87,  S21  (1990). 

I.  PULSED  FINITE  AMPLITUDE  SOUND  BEAMS 

A  time  domain  computer  algorithm  for  solving  the  KZK  (Khokhlov-Zabolotskaya- 
Kuznetsov)  nonlinear  parabolic  wave  equation  has  been  developed  for  pulsed,  ax- 
isymmetric,  finite  amplitude  sound  beams  in  thermoviscous  fluids.  This  work  has 
been  performed  by  Lee. 

A.  Background 

The  KZK  equation  accounts  consistently  for  the  combined  effects  of  nonlinearity, 
diffraction,  and  thermoviscous  absorption  on  sound  beams  radiated  by  directive 
sources.  Computer  algorithms  that  are  currently  used  to  solve  the  fully  nonlin¬ 
ear  KZK  equation  employ  the  Fourier  series  method  that  was  developed  originally 
by  Aanonsen  et  al.,1,2  and  all  computations  are  thus  performed  in  the  frequency 
domain.  The  frequency  domain  algorithms  have  been  used  with  great  success  for 
monofrequency  and  bifrequency  sources.  However,  els  either  the  source  amplitude 
or  the  number  of  source  frequencies  is  increELsed,  more  harmonics  must  be  retained 
in  the  Fourier  series,  and  the  computation  time  increases  accordingly.  Baker  and 
Humphrey,3  and  Neighbors  and  Bjprnp,4  have  recently  used  the  Aanonsen  algo¬ 
rithm  to  model  the  nonlinear  propagation  of  pulsed  sound  beams.  To  reduce  the 
number  of  harmonics  required  to  adequately  describe  a  single  pulse,  they  employ  an 
infinite  sequence  of  pulses  spaced  sufficiently  far  apart  in  time  that  adjacent  pulses 
do  not  interact  with  each  other.  Again,  however,  increasing  the  complexity  of  the 
pulse  shape  generally  requires  more  harmonics  and  therefore  greater  computation 
time,  By  performing  all  calculations  in  the  time  domain,  our  algorithm  avoids  the 
limitations  imposed  by  using  a  truncated  Fourier  series  in  frequency  domain  cal¬ 
culations.  Further  background  and  motivation  for  this  project  are  covered  in  the 
First  Annual  Summary  Report5  for  ONR  Grant  N00014-89-J-1003. 
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B.  Results 


The  numerical  solutions  of  the  KZK  equation  are  obtained  from  a  computer  al¬ 
gorithm,  developed  by  Lee,  that  performs  all  calculations  in  the  time  domain.  A 
transformation  is  applied  to  the  KZK  equation  in  order  to  facilitate  calculations 
in  the  farfield,6  rind  the  resulting  parabolic  equation  is  integrated  once  analytically 
with  respect  to  time  to  yield 


dP 

da 


1 


4(1  +  <r)2 


/( 


Vi 


p)  dr  +  + 


dr 2 


BP  dP 
(1  +  a)  dr 


(1) 


where  P,  a ,  and  r  are  dimensionless  quantities  for  the  pressure,  range  from  the 
source,  and  retarded  time,  respectively.  The  Laplacian  in  the  diffraction  term 
operates  in  the  plane  that  is  transverse  to  the  axis  of  the  beam,  .4  is  an  absorption 
coefficient,  and  B  is  a  nonlinearity  coefficient.  Finite  difference  methods  are  used 
to  integrate  the  diffraction  and  absorption  terms  over  each  incremental  step  from 
a  to  a  -f  Act.  Nonlinearity  is  included  by  changing  the  phase  r,  of  the  zth  point  on 
the  waveform  to  T{  —  BP  In  [1  +  Aa/(1  -f  <r)],  which  yields  an  exact  solution  of  the 
above  equation  when  the  diffraction  and  absorption  terms  are  ignored  ( V^P  =  0 
and  A  =  0).  Weak  shock  theory  may  also  be  included  in  the  algorithm,7  and 
focusing  may  be  taken  into  account  through  the  source  condition.8 

Numerical  results  for  the  self-demodulation  of  a  tone  burst  radiated  by  a  uni¬ 
form  piston  are  shown  in  Fig.  1.  In  the  left  column  are  shown  the  time  waveforms 
along  the  axis  of  the  beam,  and  in  the  right  column  are  shown  the  frequency 
spectra  that  are  obtained  by  taking  the  Fourier  transform  of  the  corresponding 
waveforms.  The  dimensionless  range  a  is  the  ratio  of  the  range  along  the  axis  of 
the  sound  beam  to  the  Rayleigh  distance  of  the  tone.  The  source  strength  used 
in  the  example  would  cause  the  tone  to  form  a  shock  at  a  =  1.0  in  the  absence 
of  both  diffraction  (i.e.,  for  a  plane  wave)  and  absorption.  The  absorption  length 
a"1,  where  a  is  the  attenuation  coefficient  for  the  tone,  is  a  —  0.5  in  terms  of 
the  dimensionless  range  variable.  As  the  pulse  propagates  into  the  farfield,  energy 
is  shifted  downward  toward  a  frequency  that  characterizes  the  distorted  envelope 
of  the  pulse,  in  general  agreement  with  the  pioneering  asymptotic  predictions  by 
Berktay.9  Future  work  will  involve  comparisons  of  the  results  for  self-demodulation 
with  asymptotic  predictions  not  only  by  Berktay,  but  more  recently  by  Frpysa  et 
al. 10,11  However,  the  numerical  results  are  not  subject  to  several  restrictions  that 
apply  to  the  theoretical  analyses,  e.g.,  farfield  and  weak  nonlinearity.  Additional 
investigations  will  include  strong  nonlinearity  relative  to  absorption,  and  frequency 
modulation  of  the  tone. 
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Figure  1:  Numerical  simulation  of  self-demodulation  of  a  tone  burst. 


II.  NONLINEAR  ACOUSTIC  WAVE  PROPAGATION  IN 
A  LIQUID  LAYER 

This  theoretical  investigation  was  performed  in  collaboration  with  Zabolotskaya. 

A.  Background 

In  contrast  to  the  propagation  of  finite  amplitude  sound  in  unbounded  homoge¬ 
neous  media,  the  nonlinear  propagation  of  sound  in  waveguides  and  layered  media 
is  strongly  affected  by  dispersion.  The  dispersion  law  is  determined  by  the  bound¬ 
ary  conditions  at  the  interfaces  and  the  inhomogeneity  of  the  medium,  and  it  can 
substantially  alter  the  processes  of  harmonic  generation  and  waveform  distortion. 
Natural  acoustic  waveguides  include12,13  shallow  water  channels,  isovelocity  chan¬ 
nels  in  deep  ocean,  layers  in  the  earth,  and  atmospheric  waveguides. 

In  the  present  investigation,  the  nonlinear  propagation  of  two  dimensional  sound 
waves  in  a  lossless  and  homogeneous  liquid  layer  are  considered.  The  layer  is 
bounded  below  by  a  rigid  surface  and  above  by  a  free  surface,  and  it  thus  approx¬ 
imates  the  conditions  in  a  shallow  water  channel.  The  free  surface  prevents  the 
mode  coupling  that  can  lead  to  shock  formation  in  rectangular  ducts  with  rigid 
walls.  Investigations  of  wave  distortion  and  shock  formation  in  ducts  with  rigid 
walls  are  the  subjects  of  many  papers  that  have  been  published  during  the  past 
two  decades  (see,  e.g.,  the  list  of  articles  in  Ref.  14). 

Previous  papers15-19  on  finite  amplitude  sound  in  shallow  water  channels  are 
concerned  primarily  with  parametric  (difference  frequency)  amplification  and  ra¬ 
diation  due  to  the  interaction  of  two  tones.  In  Refs.  15  and  16  the  primary  waves 
are  assumed  to  propagate  in  single  modes,  and  a  main  objective  of  the  investiga¬ 
tions  is  to  identify  resonance  conditions  for  efficient  parametric  amplification.  In 
Refs.  17-19  it  is  assumed  that  the  primary  waves  form  highly  collimated  beams 
which  are  absorbed  before  diffraction  effects  become  significant. 

B.  Results 

The  following  results  have  been  excerpted  from  Ref.  24. 

First,  a  solution  is  derived  for  second  harmonic  generation  due  to  the  propaga¬ 
tion  of  a  tone  in  a  single  mode.  The  linear  solution  which  satisfies  the  boundary 
conditions  for  a  free  surface  at  x  =  0  and  a  rigid  surface  at  x  =  h  may  be  written 
in  the  form 

pi  =  p0sin(7mx)e2!wl_fc,z) ,  (2) 

where  po  is  the  characteristic  source  pressure,  -ym  =  (2m  —  l)n/2h  defines  the 
mode  shape,  kz  =  (u  /co)\J\  —  (u>m/u;)2  is  the  propagation  wavenumber,  Co  is  the 
small  signal  sound  speed,  and  um  =  7 mc0  is  the  cutoff  frequency  for  the  mth  mode 
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(m  =  1,2,...).  We  consider  separately  the  particular  and  homogeneous  solutions 
for  the  second  harmonic  component,  and  we  write 

Pi  =  (P2)p  +  (P2)h-  (3) 


The  particular  solution,  with  pi  given  by  Eq.  (2),  is 


2  2 

(  \  - 

(P2 )?  -  _  r2,  ,2 

PO^O^rr 


Po^2  [/cq  D  \  .  2 

2  2~  (  ~7  9  J  )  Sm  (Tm2") 

[Vi  2.4/ 

y-(l  +  ^)(x  “  h)sm(2-ymx) 


where  po  is  the  static  density,  D/A  is  a  ratio  associated  with  the  quadratic  non 
linearity  in  the  isentropic  equation  of  state  that  relates  the  pressure  and  density, 
and  cz  =  u>/kz  is  the  phase  speed  in  the  2  direction.  The  mode  shapes  for  the 
forced  second  harmonic  wave  (p2)P  are  shown  in  Fig.  2  (solid  lines).  Also  shown 
are  the  corresponding  mode  shapes  for  the  fundamental  wave  pi  (dashed  lines), 
which  propagates  in  water  ( D/A  =  5),  in  the  mth  mode,  at  twice  the  cutoff  fre¬ 
quency  ( u  =  2u?m).  For  m  >  1,  second  harmonic  generation  is  seen  to  be  most 
efficient  near  the  free  surface. 

The  homogeneous  solution  is  of  the  general  form 


(Pa)h  =  £  6nsin(7nx)ej2M~Kn2) 

n=l 


(5) 


where  k„  =  (u/cq)^!  -  (u;„/2w)2  is  the  propagation  wavenumber  for  the  /ith  mode. 
The  summation  includes  both  propagating  (u;„  <  2u> )  and  evanescent  (wn  >  2~ ) 
waves.  The  propagating  waves  have  different  phase  speeds  (u;/Kn,  where  n  = 
1,2,...),  not  one  of  which  is  equal  to  the  phase  speed  of  the  forced  wave  (c,  =  u >/k. ). 
The  constants  bn  are  determined  by  the  source  condition,  which  is  assumed  to  be 
p2  =  0  at  2  =  0. 

Shown  in  Fig.  3  are  propagation  curves  and  mode  shapes  for  |p2|  obtained  from 
Eqs.  (3)-(5)  with  m  =  2,  u  =  2wm,  and  D/A  =  5.  Range  is  measured  in  terms  of  the 


Figure  2:  Mode  shapes  for  p(  (dashed  lines)  and  (p2)p  (solid  lines). 
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wavelength  A  =  2nc0/u  at  the  fundamental  frequency.  The  homogeneous  solution 
(p2  )h  accounts  for  unforced  (i.e.,  free)  sec  >nd  harmonic  waves  that  propagate  in  the 
six  lowest  modes  (1  <  n  <  6)  and  beat  with  the  forced  wave  ( p2  )p .  The  evanescent 
waves  described  by  the  homogeneous  solution  (n  >  7)  are  required  to  satisfy  the 
source  condition  p2  —  0  at  2  =  0.  Dispersion  causes  the  mode  shapes  to  vary  with 
range. 

We  now  consider  the  asymptotic  properties  of  a  slowly  varying  envelope  that 
modulates  a  high  frequency  tone.  It  is  assumed  that  the  fundamental  wave  may 
be  described  by 

Pi  =a(€r,e2z) sin(7mx)e-'M_fcl*) ,  (G) 

where  a(c-r,  e2z)  is  a  slowly  varying  envelope  with  magnitude  of  order  e.  and  r  = 
t  —  z/cg  is  a  retarded  time  based  on  the  group  velocity  cg  =  du;  fdkz.  The  envelope 
varies  more  rapidly  with  time  than  with  range,  as  indicated  by  the  explicit  depen¬ 
dence  on  the  combinations  er  and  e2z.  The  slow  scales  T]  =  er  and  Zj  =  e2z  may 
be  introduced  to  yield,  at  order  e3,  the  following  nonlinear  Schrodinger  equation 
for  the  pulse  envelope  function: 


k"  d2a  2 

+  Ta^  +  AlG|  a  =  0’ 


(7) 


where  kz  =  <Pkz/d a;2,  and  A  is  determined  by  the  nonlinear  interaction  of  the 
primary  and  second  harmonic  waves.  Solutions  of  Eq.  (7)  are  well  knowm25,26  and 
are  classified  differently  depending  on  the  sign  of  A/A:”.  With  A / k'^  <  0,  the 
modulation  is  said  to  be  stable,  and  the  amplitude  of  the  pulse  decays  according  to 
|a|  oc  1  !\fz.  With  A./kt  >0,  the  modulation  is  unstable,  and  the  pulse  may  evolve 
into  a  series  of  envelope  solitons.  The  solution  for  the  amplitude  of  a  single  soliton 
is  |a|  =  Posech(r/T),  where  Po  is  a  characteristic  amplitude,  and  T  =  P0-1  \J k"z  / A 
is  the  corresponding  period.  For  all  typical  liquids  we  found  A/A-”  <  0,  and  pulse 
propagation  is  therefore  stable.  For  the  special  case  where  |a|  is  finite  at  r  =  ±oo, 
kink  solitons  described  by  |aj  =  Pq  tanh( |t|/T)  can  exist  with  A / k,  <  0. 
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III.  NONLINEAR  EFFECTS  IN  ASYMMETRIC  CYLIN¬ 
DRICAL  SOUND  BEAMS 

This  project  was  completed  by  Kim  in  August  1990,  at  which  time  he  completed 
his  M.S.  thesis27  on  the  results  from  this  work.  The  following  background  and 
results  have  been  excerpted  from  his  M.S.  thesis. 

A.  Background 

Much  of  the  work  in  the  field  of  acoustics  has  concentrated  on  sound  sources  that 
generate  axisymmetric  sound  fields.  However,  problems  that  involve  nonaxisym- 
metric  sound  fields  include  propagation  and  reflection  of  sound  in  shallow  water, 
caustic  formation  in  inhomogeneous  media,  scattering  of  sound  by  sound,  and  rect¬ 
angular  sourc  used  in  medicine  and  in  underwater  applications. 

A  goal  of  his  investigation  was  to  develop  a  relatively  simple  numerical  model 
of  nonlinea’  effects  in  asymmetric  sound  beams.  Until  recently,  practical  con¬ 
straints  c  i  computation  time  have  limited  numerical  analyses  of  diffracting  finite 
amplitude  sound  beams  to  axisymmetric  source  conditions.  Specifically,  an  ax¬ 
isymmetric  sound  beam  may  be  described  with  only  two  spatial  coordinates,  i.e., 
the  cylindrical  coordinates  (r,  z),  where  2  measures  range  along  the  axis  of  the 
beam,  and  r  measures  distance  from  the  z  axis.  In  contrast,  a  nonaxisymmetric 
beam  must  be  described  with  three  spatial  coordinates.  The  computer  program 
that  forms  the  basis  of  the  numerical  analysis  is  an  extension  of  the  Aanonsen  fre¬ 
quency  domain  algorithm  (see  Section  I)  due  to  Bemtsen  and  Vefring.28  The  pro¬ 
gram  solves  the  KZK  equation  [Eq.  (1)  in  Section  I]  for  axisymmetric  sound  beams, 
and  in  principle,  the  algorithm  is  easily  modified  for  applications  to  nonaxisym¬ 
metric  sound  beams.  However,  the  computation  time  required  for  nonaxisymmetric 
beams  wouiu  be  roughly  the  square  of  that  required  for  axisymmet  ric  beams.  Work 
is  currently  underway  at  the  University  of  Bergen,  Norway,  to  reduce  the  computa¬ 
tion  time  required  to  use  the  program  for  general  nonaxisymmetric  fields.  Several 
examples  of  nonaxisymmetric  sound  fields  include  parametric  transmitting29'30  and 
receiving31,32  arrays,  oblique  transmission  of  a  beam  through  an  interface  33  scat¬ 
tering  of  sound  by  sound,34  35  and  wave  propagation  in  inhomogeneous  fluids.36 

The  numerical  computation  time  in  analyses  of  nonaxisymmetric  sound  beams 
can  be  reduced  considerably  by  reducing  the  dimension  of  the  computation  field.  To 
this  end,  a  two-dimensional  sound  field  is  used  to  study  nonaxisymmetric  nonlinear 
wave  propagation.  In  spite  of  the  physical  shortcomings  of  a  two-dimensional  sound 
field  in  relation  to  a  three-dimensional  sound  field  (i.e.,  cylindrical  as  opposed  to 
spherical  spreading),  the  two-dimensional  field  is  nevertheless  of  practical  interest. 
For  example,  cylindrical  spreading  is  often  an  appropriate  model  for  waveguides, 
ocean  acoustics,  and  the  nearfields  of  cylindrical  sources.  Moreover,  recent  studies 
of  caustic  formation  by  Marston, 37,38  and  of  underwater  propagation  of  pulses  by 
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McDonald  and  Kuperman,39  are  performed  for  two-dimensional  sound  fields.  In 
cases  where  cylindrical  spreading  is  an  inappropriate  approximation  of  a  given  finite 
amplitude  sound  field,  the  analysis  in  two-dimensions  should  still  provide  physical 
insight  into  phenomena  which  occur  in  the  corresponding  spherically  spreading 
field. 


B.  Results 

Closed  form  analytical  solutions  were  derived  for  linear  radiation  from  an  amplitude' 
shaded  piston  (two  dimensional),  and  for  second  harmonic  generation  in  a  Gaussian 
beam.  The  analytical  solutions,  which  were  derived  from  the  KZIv  equation,  were 
used  to  verify  the  accuracy  of  the  numerical  solutions,  and  to  gain  insight  into  the 
asymptotic  properties  of  cylindrically  spreading  sound  beams.  Two  problems  were 
studied  in  detail  with  the  computer  program,  high  intensity  radiation  from  ampli¬ 
tude  shaded  piston  sources,  and  scattering  of  sound  by  sound  from  two  Gaussian 
beams. 

Graphical  results  from  the  investigation  of  amplitude  shaded  pistons  were  pre¬ 
sented  in  the  First  Annual  Summary  Report.5  The  results  show  that  nonlinear 
effects  are  altered  by  asymmetric  source  conditions.  Namely,  the  definition  of  the 
sidelobe  structure  in  the  higher  harmonic  beam  patterns,  including  the  nonlinear 
nearfield  effects  known  as  fingers,6  are  very  sensitive  to  source  asymmetry.  We  also 
observed,  in  the  results  from  both  nonlinear  and  linear  theories,  that  the  degree  of 
asymmetry  in  the  beam  decreases  with  range  from  the  source. 

The  fundamental  aspects  of  scattering  of  sound  by  sound  were  found  to  be  very 
similar  in  both  spherically  spreading  and  cylindrically  spreading  primary  beams. 
This  conclusion  was  reached  by  comparing  quasilinear  numerical  results  for  scatter¬ 
ing  of  cylindrical  Gaussian  beams  with  quasilinear  analytical  results  for  scattering 
of  spherical  Gaussian  beams.35  Specifically,  the  angle  at  which  scattered  sound  is 
most  prominent  is  the  same  in  the  two  cases,  and  the  dependence  of  the  relative 
amplitude  of  the  scattered  sound  on  various  source  parameters  (e.g.,  beam  intersec¬ 
tion  angle,  frequency  ratio,  etc.)  also  follows  the  same  general  trends  in  each  case. 
Fully  nonlinear  scattering  of  sound  by  sound  was  also  investigated  numerically. 
We  observed  that  noncollinear  interaction  of  sound  beams  may  produce  highly  os¬ 
cillatory  beam  patterns  due  to  interference  between  harmonics  that  propagate 
different  directions.  Higher  order  scattering  that  is  not  predicted  by  quasilinear 
theory  was  also  observed.  In  special  cases,  however,  quasilinear  theory  was  found 
to  be  instrumental  in  predicting  the  direction  of  the  higher  order  scattering. 
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IV.  EFFECTS  OF  ABSORPTION  ON  THE  INTERAC¬ 
TION  OF  SOUND  BEAMS 

This  work  was  performed  by  Darvennes,  Naze  Tj0tta,  and  Tjotta  as  an  extension 
of  earlier  research  on  scattering  of  sound  by  sound,  as  discussed  in  the  First  Annual 
Summary  Report.5  This  project  is  now  completed. 

A.  Background 

In  the  previous  annual  summary  report,  the  effects  of  absorption  on  the  nonlinear 
interaction  of  two  Gaussian  beams  were  discussed.  The  results  of  that  work  were 
reported  in  Ref.  40.  During  the  past  year,  the  numerical  results  were  extended  to 
include  second  harmonic  generation  in  a  sound  beam  radiated  by  a  circular  piston. 

B.  Results 

The  following  summary  is  excerpted  from  a  paper41  that  is  currently  in  press. 

Asymptotic  formulas  were  derived  for  the  sum  and  difference  frequency  sound 
that  is  generated  by  the  interaction  of  two  harmonic  sound  beams  in  a  dissipa¬ 
tive  fluid.  The  results  are  valid  for  arbitrary  source  configurations,  within  the 
limits  of  the  parabolic  approximation.  At  range  2,  different  asymptotic  results 
apply  to  the  cases  of  moderate  absorption  [z  —  0(\L±\)  »  R±\  and  strong  absorp¬ 
tion  (2  »  R±  »  \L±\),  where  L±  =  (oq  +  a2  —  is  an  absorption  length, 

R±  =  2klk2a2 /k±  is  a  diffraction  length,  k  and  a  and  the  wavenumbers  and  ab¬ 
sorption  coefficients,  respectively,  at  the  primary  frequencies  uq  and  u>2  and  the 
sum  and  difference  frequencies  ui±  —  uq  ±  u>2,  and  a  is  the  characteristic  radius  of 
the  source.  It  was  found  that  the  relative  amplitudes  of  the  nonlinearly  generated 
pumped  and  scattered  waves  (see  Refs.  34  and  35  for  precise  definitions  of  the 
pumped  and  scattered  waves)  depend  not  only  on  the  nature  and  geometry  of  the 
sources  (e.g.,  ka  of  each  source,  on-source  pressure  distributions,  source  separation, 
and  beam  interaction  angle),  but  also  on  the  dissipative  properties  of  the  medium 
through  the  sign  of  L±.  In  the  farfield,  the  scattered  sound  dominates  the  pumped 
sound  at  combination  frequencies  for  which  L±  >  0,  whereas  the  reverse  is  true  at 
combination  frequencies  for  which  L±  <  0.  Numerical  results  have  been  obtained 
for  the  case  of  two  Gaussian  beams  that  cross  at  an  angle,  and  for  the  second  har¬ 
monic  component  generated  by  a  circular  uniform  piston  source.  Various  frequency 
dependencies  were  considered  for  the  absorption  law:  quadratic,  linear,  and  square 
root. 
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V.  PARAMETRIC  RECEIVING  ARRAYS 


This  project,  which  is  now  completed,  was  performed  in  collaboration  with  Dar- 
vennes  and  Naze  Tjptta. 

A.  Background 

The  motivation  for  this  research  is  the  potential  use  of  a  parametric  receiving  ar¬ 
ray  for  measuring  freefield  source  directivities  in  reverberant  environments.  Back¬ 
ground  on  this  project  may  be  found  in  Ref.  31. 

B.  Results 

New  results  during  the  past  year  have  included  the  effects  of  absorption  on  para¬ 
metric  reception.  The  following  summary  is  excerpted  from  a  paper42  that  was 
presented  at  the  12th  International  Symposium  on  Nonlinear  Acoustics  in  Au¬ 
gust  1990. 

Beam  patterns  are  obtained  with  a  parametric  receiving  array  (parray)  by  ro¬ 
tating  a  high  frequency  (uq)  pump  beam  through  the  entire  low  frequency  (u.-2 ) 
wave  field,  with  the  hydrophone  always  located  on  the  axis  of  the  pump  beam.  For 
our  application,  the  pump  is  located  at  the  center  of  a  low  frequency  source,  and 
the  sum  and  difference  frequency  pressure  p±(0,z )  is  recorded  at  the  hydrophone. 
When  there  is  no  absorption  and  no  reflector,  the  pressure  p±(8,z)  exhibits  an  an¬ 
gular  dependence  that  converges  as  z  — >  oo  to  the  product  of  the  freefield  directivity 
functions  of  the  primary  beams,  DX{8)D2{0).  Given  the  known  pump  directivity 
Di(0),  an  estimate  of  the  unknown  source  directivity  D2(8)  may  thus  be  obtained 
from  the  normalized  pressure  P±(0,z)  =  p±(0,  z)/Di(8).  In  the  absence  of  absorp¬ 
tion,  however,  P±{8,z )  is  a  good  approximation  of  the  farfield  directivity  function 
D2(6)  only  when  In (z/zj)  >  1  (where  j  —  1,2,  and  Zj  is  the  Rayleigh  distance  of 
the  jth  primary  beam,  i.e.,  of  either  the  low  frequency  source  or  the  pump).  In 
the  presence  of  a  reflector  and/or  of  dissipation,  the  relation  between  P±(0,z)  and 
D2(0 )  may  be  more  complex. 

As  discussed  in  the  previous  annual  summary  report,5  the  parray  is  capable  of 
discriminating  against  multipath  interference  due  to  reflections  from  an  adjacent 
surface.  However,  the  efficiency  of  the  discrimination  depends  on  the  distance 
between  the  pump  and  the  hydrophone.  As  the  distance  is  increased,  the  multipath 
components  become  increasingly  significant.  When  the  hydrophone  is  close  to 
the  source,  the  interference  less  significant,  but  P+(0,  z)  slightly  overestimates  the 
beamwidth  of  D2(0).  The  overestimation  of  beamwidth  due  to  finite  hydrophone 
ranges  is  discussed  in  detail  in  Ref.  6. 

In  Fig.  4  we  show  the  effect  of  thermoviscous  absorption  on  the  normalized 
sum  and  difference  frequency  pressures  P±(8,z).  The  absorption  increases  from 
left  to  right  in  the  figures,  with  Fig.  4(a)  showing  the  case  of  no  absorption.  The 
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Figure  4:  Effect  of  absorption  on  parametric  reception  near  a  reflecting  surface. 


sum  frequency  beam  pattern  P+(#,2)  is  the  solid  line,  the  difference  frequency 
beam  pattern  P_(#,2)  is  the  dot-dash  line,  and  the  low  frequency  beam  pattern 
(that  which  is  to  be  measured  with  the  parray)  is  the  dashed  line.  The  range  is 
z  —  5022  and  the  source  depth  is  given  by  h  =  2e2,  where  e2  is  the  source  radius, 
and  k\ / k2  =  100,  £2/e\  =  10,  k2e2  =  30.  As  the  absorption  parameter  |c*t|~2  is 
increased  from  0  to  0.1,  the  sum  frequency  directivity  P+(9,z)  tends  toward  the 
low  frequency  beam  pattern.  In  contrast,  for  large  values  of  absorption,  P_(0,2) 
tends  toward  a  different  directivity  function  which  depends  primarily  on  the  ratio 
k-/otT.  The  sum  frequency  signal,  rather  than  the  difference  frequency  signal, 
therefore  provides  more  information  about  the  low  frequency  source  distribution 
when  absorption  is  significant.  This  conclusion  also  follows  from  the  more  general 
analysis  of  absorption  in  Section  IV  of  this  report. 
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